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Schwinger bosons allow for an advantageous representation of quantum double- exchange. We review 
this subject, comment on previous results, and address the transition to the semiclassical limit. We 
derive an effective fermionic Hamiltonian for the spin-dependent hopping of holes interacting with 
a background of local spins, which is used in a related publication within a two-phase description 
of colossal magnetoresistant manganites. 

PACS numbers: 75.10.-b General theory and models of magnetic ordering, 75.30.Et Exchange and 
superexchange interactions, 75.30.Vn Colossal magnetoresistance 
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I. INTRODUCTION 

Introduced by Zener p| in the early 1950s the notion 
of double exchange together with mixed- valency mangan- 
ites Ri^ x AJAvlO z (where R = La, Pr, Nd and A = Sr, 
Ca, Ba, Pb) attracted renewed attention when a colos- 
sal magnetoresistive effect was discovered in these com- 
pounds some years ago 0. The magnetic and electronic 
properties of manganese oxides, to some extend, are be- 
lieved to arise from the large Coulomb and Hund's rule 
interaction of the manganese d shell electrons. Due to 
the almost octahedral coordination within the perovskite 
structure the d levels split into two subbands labeled ac- 
cording to their octahedral symmetry, e g and t2 g . In the 
case of zero doping (x = 0) there are four electrons per 
Mn site which fill up the three t^g levels and one e g level, 
and by Hund's rule, form aS = 2 spin state. Doping will 
remove the electron from the e g level, and by hopping via 
bridging oxygen sites these holes acquire mobility. How- 
ever, this hopping acts in a background of local spins 
5 = 3/2 formed by the ti g electrons and its amplitude 
depends on the overlap of the spin states at neighbouring 
sites (or, in a classical language, on their relative angle), 
it is largest if the total bond spin is maximal and vice 
versa ||. 

Another ingredient, that is assumed to significantly in- 
fluence the physical properties of manganites, is electron- 
lattice interaction. Namely the two e g orbitals, which are 
degenerate in a perfect cubic environment, will couple to 
lattice vibrations of the same symmetry, giving rise to 
a Jahn- Teller effect and polaronic behaviour in some re- 
gions of the phase diagram. It is this close interplay of 
three different subsystems (electrons in degenerate or- 
bitals, background of localized spins, and lattice vibra- 
tions) that makes the physics of manganites both, rich 
and complicated. 

In the present work we concentrate on the double ex- 
change (DE) part of the interactions and consider dif- 
ferent possibilities for an approximate treatment of the 
exact DE Hamiltonian on a lattice in terms of effective 
electronic one- or two-band models. These can be used in 



a more elaborate modelling of the real materials (see our 
forthcoming work Q). It turns out that quantum double 
exchange on a lattice is most suitably derived and de- 
scribed with the help of Schwinger bosons. We therefore 
include a detailed and pedagogic derivation of the quan- 
tum DE Hamiltonian using Schwinger bosons. Although 
this approach has been used before (5||(|, we feel a com- 
prehensive presentation of the subject is still missing. In 
two appendices we reexamine the derivation for two sites, 
and consider the semiclassical limit (S — + oo). In addi- 
tion, by means of numerical experiments, we illustrate 
how this limit evolves from the quantum case. 



II. SCHWINGER BOSON REPRESENTATION OF 
DOUBLE EXCHANGE 

To derive the quantum DE Hamiltonian on a lattice, 
as a starting point we take the Kondo lattice model in- 
cluding on-site Coulomb repulsion, 
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where summation is over nearest neighbour bonds (ij) 
or sites i, respectively. For clarity and since it can be 
included easily in the final result, here we have neglected 
the orbital degeneracy of the e g electrons. That is, c^J 
denote electrons in a single band, which interact with 
some localized spin Sj via the Hund's coupling Jh- In 
the real materials this localized spin corresponds to the 
remaining t2 g electrons, that tend to form a high spin 
state with an electron in the e g shell. 

In the manganites the situation is such, that U 3> 
Jh > t (cf. Refs. pj7| ). Hence, we first take the limit 
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U — > oo, resulting in 



with the projectors 



H = -t 
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with restricted fermions Ci a = Ci a (\ — n;_ CT ). Next, fol- 
lowing Kubo and Ohata || , the exchange term in Eq. 
is solved while the hopping term is considered as a small 
perturbation. For positive Jjj the ground-state of the 
exchange term is a free spin 5, if there is no electron at 
site i, or a coupled spin 5 = 5+1/2 otherwise (note, that 
we use 5 for the length of the localized spin Si formed 
by t2 g electrons). To describe the effective hopping we 
therefore need a projection operator, which restores these 
conditions, 



(s t(W ) + (g + i)a g 
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Then the DE-Hamiltonian (Eq. (2.3) in Ref. g) is given 
in terms of spin and restricted fermion operators by 



C = -*E [4(^ + p+u,^, + h.c. 



(4) 
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However, this expression turns out to be unwieldy for an- 
alytic as well as numeric calculations. Although in prin- 
ciple the electronic spin is absorbed into the total spin 
at each site, the spin index a is still present in Eq. (Q). 
Here the advantages of Schwinger bosons come into play, 
namely the possibility to describe spins of arbitrary am- 
plitude with the same set of boson operators et^ and bi , 



5+ = a\b l , S i = b\ai , 
Sf = (oJo < -6t6.)/2, 
\S i \ = iata i +btb i )/2. 
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Using these operators, we can rewrite the projection op- 
erators Pi , 
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where we can keep 5 in the denominator, because it is 
conserved. The last matrix can be decomposed easily, 



which leads to 



1 



25 + 1 



(9) 



prDE _ 

"25 



■t 



H.c. 



(10) 



R 



V2S~- 



(11) 



If we restrict the Hilbert space of the problem by fixing 
the spin length at each site to the value 5 + 77,^/2 (where 
rii = rii-\ + ^u)j we can simply replace the projectors 
R~l by spinless fermion operators c,. This can be seen, 
by taking a closer look on the operators Rf . Let us 
assume, site i is occupied by an electron, and therefore 
the total spin at the site is 5 + 1/2. The corresponding 
state \n e ] 5, m) can be written as 



\l;S+km) 
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Using the representation of an 5 Z eigenstate in terms of 
Schwinger bosons, 
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we find, that applying the operator R + the state |1; 5 + 
1/2, m) in Eq. ( p"2| ) is transformed into the corresponding 
Schwinger boson representation of the coupled spin 5 + 
1/2, while the electron is annihilated 



B+\l\S+\,m) = |0;5+ \,m) 



(14) 



Backwards, the operator (R + Y creates the decomposi- 
tion of the coupled spin into electronic and localized spin, 
i.e., it produces appropriate Clebsch-Gordan coefficients. 

Now it is straightforward to omit the electronic spin 
index, using only spinless fermions, , and Schwinger 
bosons to describe double-exchange. The corresponding 
Hamiltonian is given by 
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where, for every site i, the Hilbert space is constrained 
to 



a\a t + b\b t = 25 + cjc, : 



(16) 



In the case of low doping usually it is more appropriate 
and natural to consider holes instead of electrons. Here 
"hole" denotes a fermion and a spin 5 moving together 
in a background of spins 5 = 5+1/2. For the trans- 
formation of the electronic Hamiltonian, Eq. (|l0|), the 
operator R + needs to be replaced by a suitable counter- 
part involving restricted hole operators hi a . In analogy 
to Eqs. dill) to dll), the state 
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is transformed into its corresponding Schwinger boson 
representation by the operator 



htfbi — hiicti 



(18) 
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Hence, in hole representation, the DE Hamiltonian reads 
t 
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£ [(RT)H*i<4 + t>$)Rj + H.c] . (19) 



The Hamiltonian for spinless fermions, Eq. (|15|), changes 
only little, becoming 



with the constraint 



a} ai + b}bt=2S-h]h 



(20) 
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III. EFFECTIVE TRANSPORT HAMILTONIAN 

To obtain an effective Hamiltonian for the spin- 
dependent hole-hopping, the spin part of the DE interac- 
tion is considered within mean field approximation. How- 
ever, there are two representations of the DE Hamilto- 
nian to start from: Eqs. (19) and (|2(i|). The resulting 
effective Hamiltonians describe carriers with or without 
spin, respectively. Given Eq. (|20|), a mean hopping of 
the spinless carriers, t^ h \ is obtained by considering each 
bond (ij) separately, and taking the expectation value of 
the DE term for all values of the total bond spin St and 
S T in an effective ordering field A = PgfisH^ s . This re- 
produces the result of Kubo and Ohata || . Having spin 
S at site i and spin S at site j, the coupled state St with 
maximal S T is given by 

\St, St)(§s) = 

= C (6to] - 6t a |)2S+i-ST (a t ) s T +i (a t ) s T -i| 0) _ (22) 
where C is a normalization factor, 
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Applying (aja i + b] j b i ) we arrive at 



( a ] a i + b ] b i)\ S T,S T )(§S) = 

= C(S T + !)(&tat - b]4f s +i-^(al)^-Ha}f T ^\0) 
= (St + ^\St,S t ) ( ss)- (24) 

Hence, we rederived the effective matrix element for a 
single bond (cf. Appendix |X] or Ref . || ) , 
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(b) _ 5j_+_l/2 
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(25) 



which is averaged over all values and directions of St, 
m= 7§ [SX]t, (26) 

where [pi 



7s [SX] = 



2S-1/2 
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EOT 
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2S e 
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S r =l/2 2^iM=-Si 
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_coth(2^hiA) [coth(S*A) - icoth(|)] 



The effective Hamiltonian describing spinless fermionic 
holes in an averaged background of ordered spins reads 



rreff.1 _ r(b) \p 
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In Appendix [b| we compare the classical limit of both, 
this Hamiltonian and the exact expression. At least the 
bandwidth turns out to be represented very well. 

Another possible way to obtain an effective Hamilto- 
nian is based on the picture of itinerant carriers of spin | 
moving in the background of localized spins, the correla- 
tions of which change on a large time scale compared with 
the hopping frequency. Then an effective h opp ing Hamil- 
tonian is obtained averaging H^ c , Eq. (|19|), over free 
spins S in a homogeneous field A. According to Eqs. @- 
(§) and the fact that (S ,± ) = and (S z ) = SB S [SX] 
(where -Bg[z] denotes the Brillouin function), only two 
terms contribute. The resulting Hamiltonian involves 
two effective hopping matrix elements, one for each spin 
channel, 



h: 



cff,II 
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H.c. 



with 



(29) 

(30) 
(31) 



In a fully polarized background (A — > oo) only holes with 
anti-parallel spin can hop (ti — > jf+T i while holes with 
parallel spin are blocked (tf — » 0). In general, the situ- 
ation is complicated by the fact that Eq. ( p9| ) involves 
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restricted fermion operators (Hubbard operators) forbid- 
ding an exact solution of the model. However, to a good 
approximation, in the polarized phase the up-band can 
be neglected, while the down-band is taken into account 
using unrestricted operators h^V . On the other hand, in 
a disordered phase (A — > 0) both bands are equivalent 
(tf = ti — > 2S+1 ma k m g an approximate treatment 
less evident. 



IV. CONCLUSIONS 

In the present work we review the subject of double ex- 
change and derive an effective Hamiltonian for the spin- 
dependent hopping of holes in an averaged background 
of local spins. This is used in a subsequent publication 
within a two-phase scenario for the description of colossal 
magnetoresistant manganites, alternatively to the effec- 
tive Hamiltonian for spinless carriers derived in Ref. 
Besides, we illustrate that on the level of quantum spins 
all features of double exchange - namely its derivation, 
limiting cases, as well as approximations - are probably 
most clearly represented in terms of Schwinger bosons. 
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APPENDIX A: REEXAMINATION OF THE 
2-SITE PROBLEM 



As noted above, the DE matrix element, Eq. (f25[), was 
first derived by Anderson and Hasegawa @] by consider- 
ing a system of two spins S at neighbouring sites and a 
mobile electron whose spin a is coupled to the local spin 
at the same site. Recently the problem was reexamined 
by Miiller-Hartmann and Dagotto || , who found a "non- 
trivial" phase factor, which the authors interpreted in 
terms of a Berry phase within the limit S — > oo. We argue 
that the phase factor of their quantum DE Hamiltonian 
merely compensates for the phase introduced by the per- 
mutation of neighbouring spin states, whereas indepen- 
dently the Berry phase evolves in the classical (S — > oo) 
limit of the correct DE Hamiltonian, Eq. Jl5|). More- 
over, on a lattice double exchange can not be described 
in terms of spin and permutation operators, which do not 
take into account fcrmionic commutation relations. 

To be specific we briefly summarize the two site prob- 
lem. Assuming the electron initially to be at site 1 we 
start with coupling the local spin Si and the electron spin 
<t, to give an on-site spin Si = Si +cr. We note, that the 
construction of the corresponding wave function is not 
unique, as is known from textbooks on spin algebra [f[o|| . 
In more detail 
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where we give the order of the constituting spin states 
as a subscript and denote Clebsch-Gordan coefficients by 
(•• | •). Adding this spin to the unchanged core spin at 
site 2, we arrive at an initial state of total spin St 



\<p b 



|S , T m T ) ( ( cr 5 1 ) 52 

E 



{ B1 


s 2 


St \ 




m 2 


m T J 



\Simi) 



(<rSi) 



\S 2 m 2 ) 



(A3) 
(A4) 



Now, if the electron moves to site 2, we can couple <x 
and S2 to give §2, which together with Si yields another 
state of total spin St- However, as can be seen above, 
there are different ways to connect the three spins. One 
possible final state is 



\S T mT)(Si(aS2)) 

E 



(ft 


s 2 


S T \ 


\ mi 


m 2 





|SiTOi)|S 2 TO 2 ) 



(<rS 2 ) 



(A5) 
(A6) 



(b) 



The corresponding effective hopping matrix element t A 
is proportional to the overlap (ip^\ip in ), which can be 
evaluated by rewriting the multiple sum over the product 
of four Clebsch-Gordan coefficients with a 6-j-symbol, 



f_-j\Si+S 3 +S T 



(A7) 



(25i + 1)(2S 2 + 1) 



a Si Si 
St S 2 S 2 



We obtain the value of the matrix relevant for strong 
Hund's coupling by setting |Sj| = Si + 1/2 = S and 
|S<| = S = S- 1/2: 
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2S 
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i.e., there is no Sr-dependent phase factor. However, 
Miiller-Hartmann and Dagotto derive an effective double- 
exchange Hamiltonian, where the hopping is expressed as 
the permutation of the spin S "particle" at site 1 and the 
spin (S — 1/2) "hole" at site 2. Therefore the final state 
they have to consider for the matrix element is 



\(fi) = |Stto t ) ((ctS2)Si) 
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Obviously the permutation of the indices in the Clebsch- 
Gordan coefficients yields a different phase factor in the 
overlap 



(25i + l)(2S 2 + i; 



(All) 

a Si Si 1 
St S 2 S2 J 
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and consequently the authors obtain a Sr-dependent 
phase factor for ta , 



(A12) 

To recover the effective Hamiltonian Eq. (6) of Ref . [|| , 
H cS = -t P 12 Q s (y) , (A13) 

we simply have to express St in with the help of 
y = Si ■ S 2 /(S(S — 1/2)). In a direct way we can use 
interpolation polynomials in the form of Lagrange, i.e. 



for a given spin configuration {6^, 4>k} and two electronic 
states | -0i } and \1jJ2), 

Hj) = I[Kkm(S +^,9 k ,<f> k )) (B4) 

k 

(where \rij.k) — (cj,)™^ |0) with numbers n^h € {0,1}), 
we find the average 



(hhtih = n< ni .*i( - e Ms- + h - c -] ) n 



n 2 ,k) 
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with the matrix clement 



2S-1/2 2S-1/2 

,(b) tt v-v* 



ys T - Vj 
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Qs(y)= E 4 b) II 

S T =l/2 j^ 1 / 2 

_ j(j + 1) - S(S + 1) - (S - 1/2)(S + 1/2) 
% ~ 25(5-1/2) ' ( j 

or the recursive formula given in Eq.(7) of Ref. |J. Of 
course, taking the matrix element for the construc- 
tion of Q§(y) is wrong, as this does not account for the 
phase factor due to the permutation P 12 . Concerning 
the limit of classical spins, we can see no connection 
between the above Sr-dependent phase factor and the 
Berry phase. 

The peculiarities concerning the phase factor indicate 
that the above derivation is not suitable for the general- 
ization to a lattice. Moreover, the spin-5 "particles" still 
obey fermion commutation relations, which can not be 
expressed by conventional permutation or spin operators. 
That means, the expression given in Eq.(6) of Ref. 



tj, = cos I — I cos I — l e 
13 I 2 J \ 2 



-i(0i-^-)/2 



2 \ e i(4>i-<t>j)/2 
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Hence, the classical Hamiltonian should read 



# eff = -tJ2 p ijQs(y)' 



(Ai6) 

is not the correct quantum DE-Hamiltonian on a lattice. 



APPENDIX B: CLASSICAL LIMIT OF THE DE 
MODEL 

The limit S -> 00 of H^ E , Eq. (|l|), is easily de- 
rived by taking its expectation value with spin coherent 
states jn|, 



|0(W)) = M ^ )2S 10), 



(Bl) 



where u = cos(#/2) e^ 2 and v — sin(#/2) e l <t>/ % . Using 
the properties of coherent states, 

a \(l(S t 9,4)) = V2Su |0(5- ±,0,0)) (B2) 

b |n(S,0,0)) =V2Sv \n(S- |,M)>. (B3) 
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(B7) 



which is equivalent to the results obtained in Refs. Jl2],p[. 
In our case, however, the classical limit followed from the 
quantum Hamiltonian, Eq. ([l5]), in an obvious and more 
straightforward way. Note also, that averaging the oper- 
ators (Rf Y and Rf over coherent states yields the uni- 
tary transformation onto rotated electrons d\ (compare 



Refs. |2JB) 



di = cos 



I 1 



Si ' 



(B8) 



i.e., naturally Eq. ( fl0| ) has the same classical limit. 

0.8 




-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



FIG. 1. Density of nonzero eigenvalues (dashed line) and 
running average (bold dot-dashed), calculated for 2 electrons 
on 4 sites with S = 5 (inset) and S = 10, compared to the 
classical result 5* — > 00 (bold solid). 
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To check, whether the description of double-exchange 
in terms of classical spins is appropriate, we compared the 
(canonical) density of states (DOS) for a fixed number of 
carriers on a small cluster, whic h interact with quantum 
(Eq. (|l5|)) or classical (Eq. (B7)) spins, respectively. Us- 
ing Chebychev expansion and maximum entropy meth- 
ods for the quantum case the spectra can be ob- 
tained numerically for rather large S. The classical DOS 



is found by averaging the eigenvalues of H^f SB , Eq. (B7), 
over a large number of spin configurations. 
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E 

FIG. 2. Density of nonzero eigenvalues for 4 electrons on 8 
sites with S = 3/2; same notations as in Fig. |lj 

In Fig. [l] we consider two electrons on a ring of four 
sites. Comparing the running average (bold dot-dashed) 
over the discrete spectrum (thin dashed) and the classical 
limit (bold solid), we find good convergence already for 
a moderate spin length S = 10. But even for the case 
S = 3/2, which is realized in the manganites, the classical 
description appears to be acceptable. If we consider four 
electrons on a ring of eight sites, the spectrum is much 
more dense, making similarities to be recognized easily, 
see Fig. [|. Of course, from the density of states we can 
learn nothing about correlations or other more involved 
features. Note, that in both figures we subtracted the 
peak at E = consuming a large fraction of spectral 
weight. 

Another interesting check concerns the effective hop- 
ping matrix element t^°\ Eq. (|2^). With classical spins 
and the Chebychev expansion methods mentioned above, 
one can calculate the grand-canonical DOS of the tight- 
binding model, Eq. (B7), for rather large clusters (here 
64 3 sites on a simple cubic lattice) without much effort. 
By considering a thermalized ensemble of classical spins 
in a homogeneous field we compare the resulting band- 
width with the limit S — > oo of f( h \ where the limit of 
7 i g[<S'A] is obtained omitting S in the argument and set- 



ting S — > oo in the index 
1 



7s- 



>[A] 



1 + coth(A) 



coth(A) 



(B9) 



The agreement is rather satisfactory, as can be seen in 
Fig. ||. Naturally, the precise shape of the DOS (see inset) 
will not be reproduced by the effective electronic model, 
it remains simple cubic tight-binding for all fields and 
temperatures. 




FIG. 3. Bandwidth of tight-binding electrons on a simple 
cubic lattice of size 64 3 , interacting with thermalized clas- 
sical spins in a homogeneous field A, compared to the limit 
S — > oo of the Kubo/Ohata formula; inset: Narrowing of the 
corresponding density of states with decreasing A. 
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